Defect Energy Levels in Density Functional Calculations: Alignment and Band Gap 

Problem 
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TABLE I: 


Calculated and 


experimental 


band gaps (in eV). 




Si 


SiC 


HfOa SiOa 


Semilocal 

Hybrid 

Expt. 


0.6 
1.8 
1.2 


2.2 
3.9 
3.3 


4.3 5.8 
6.7 8.3 
5.9 8.9 



Audrius Alkauskas, Peter Broqvist, and Alfredo Pasquarello 

Ecole Polytechmque Federate de Lausanne (EPFL), 

Institute of Theoretical Physics, CH-1015 Lausanne, Switzerland and 

Institut Romand de Recherche Numerique en Physique des Materiaux (IRRMA), CH-1015 Lausanne, Switzerland 

(Dated: July 29, 2008) 

For materials of varying band gap, we compare energy levels of atomically localized defects calcu- 
lated within a semilocal and a hybrid density-functional scheme. Since the latter scheme partially 
relieves the band gap problem, our study describes how calculated defect levels shift when the 
band gap approaches the experimental value. When suitably aligned, defect levels obtained from 
total-energy differences correspond closely, showing average shifts of at most 0.2 eV irrespective of 
band gap. Systematic deviations from ideal alignment increase with the extent of the defect wave 
function. A guideline for comparing calculated and experimental defect levels is provided. 
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Semilocal approximations to density functional the- 
ory, such as the local density approximation (LDA) and 
the generalized-gradient approximation, have proved ex- 
tremely valuable to investigate energetic, atomistic, and 
magnetic properties of defects in solids [Ij. However, 
these approximations have been much less successful in 
locating charge transition levels in the band gap, be- 
cause of the well known band gap problem from which 
they suffer. As a result, a direct comparison between 
calculated and experimental energy levels remains am- 
biguous. Furthermore, the determination of equilibrium 
densities of intrinsic defects and charge carriers is hin- 
dered [4]. Therefore, considerable efforts have been de- 
ployed in the study of defects to address the band-gap 
problem going beyond semilocal approximations to den- 
sity functional theory. Many-body perturbation theory 
in the GW approximation is the method of choice for 
calculating defect levels [3[, but remains computation- 
ally demanding and therefore limited to small-size sys- 
tems. Several practical routes have also been proposed, 
such as the scissor-operator scheme, the marker method 
[j|, the LDA-|-f7 method [5|, the use of adapted pseu- 
dopotentials [6|], and the application of ad- hoc extrap- 
olation schemes [7|. However, the general applicability 
of these approaches is unclear. More recently, hybrid 
density functionals have become increasingly popular for 
addressing defect energy levels Q- These functionals in- 
corporate a fraction of Hartree-Fock exchange, leading 
to higher accuracy [9| and improved band gaps [10| com- 
pared to semilocal functionals. 

In this work, we carry out a comparative study be- 
tween defect energy levels calculated with semilocal and 
hybrid density functionals to determine their shifts as 
the description of the band gap improves. We aim at 
gaining insight into how calculated and measured defect 
levels should be compared when the adopted theoretical 
scheme is subject to the band gap problem. For this pur- 
pose, we considered materials covering a large range of 



band gaps and selected defect levels spanning large por- 
tions of their band gaps. Our study shows that charge 
transition levels obtained with semilocal and hybrid den- 
sity functionals correspond closely, provided a suitable 
alignment scheme is adopted. As the band gap decreases, 
systematic deviations from ideal alignment are found to 
increase with the extent of the defect wave function. 

The semilocal density-functional calculations were per- 
formed within the generalized gradient approximation 



proposed by Perdew, Burke, and Ernzerhof (PBE) [11|. 
We used a hybrid density functional, denoted PBEO, 
which is obtained from the latter by replacing 25% of the 
PBE exchange energy by Hartree-Fock exchange [I4I . We 
used a scheme based on plane-wave basis sets and norm- 
conserving pseudopotentials. The pseudopotentials were 
generated at the semilocal level and used in all calcula- 
tions. The plane- wave basis set was defined by an energy 
cutoff of 70 Ry. The Brillouin zones of our supercells were 
sampled at the F point, but primitive cells with a con- 
verged /c-point sampling were used for the determination 
of the bulk band edges. We took care of the integrable di- 



vergence of the Hartree-Fock exchange term [13| . Struc- 
tural relaxations were carried out at the semilocal level 
14| . Our calculations were performed with the codes 
Quantum-ESPRESSO [il| and CPMD flit . 

We considered defects in four different materials of 
varying band gap: Si, SiC, Hf02, and Si02. The band 
gaps of these materials calculated at the semilocal level 
severely underestimate the experimental values (Table|l|, 
as usual for semilocal density functional schemes. As 
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FIG. 1: (Color online) Schematic illustration of the alignment 
between energy levels obtained with a semilocal and a hybrid 
density functional. The charge transition levels n and Jl are 
referred to the respective valence band maxima (VBM) and 
to a common reference level, respectively. The conduction 
band minima (CBM) are also shown. 



shown in Table [H the hybrid scheme systematically gives 
larger band gaps, generally leading to a better agreement 
with experiment. For silicon, we adopted a cubic simula- 
tion cell of 64 atoms and considered the following defects 
with their relative charge states: the Si vacancy (+2, +1, 
0, —1, —2), the Si self-interstitial (-1-1, 0), the substitu- 
tional O (0, —1), the C interstitial (-1-1, 0, —1), and the 
P/Si-vacancy complex (+1, 0, —1). For SiC, we modeled 
the AH polytype using an orthorhombic cell containing 
96 atoms. Considered defects include the C vacancy {+2, 
+1, 0, -1, -2), the Si vacancy (-^2, +1, 0, -1, -2), and 
the complex consisting of C substitutional to Si next to 
a C vacancy (-1-2, -|-1, 0, -1, -2). For Hf02, we took the 
monoclinic structure and used a supercell containing 96 
atoms. We considered the O vacancy (-1-2, —1, 0, -1-1 and 
—2) and the O interstitial (0, —1, and —2). For Si02, we 
modeled a-quartz with an orthorhombic cell containing 
72 atoms and considered the interstitial H (-1-1, 0, —1), 
the Si-Si dimer bond (4-1, 0), the puckered O vacancy 
{+1, 0), the H bridge -Si-H-Si- (-f 1, 0, -1), the substi- 
tutional N (0, -1), and the interstitial O2 (0, -1). All 
defect states studied are atomically localized. 

The formation energy of a defect in its charge state q 
can be expressed in terms of the electron chemical po- 
tential /i referred to the valence band maximum Ev [I]'- 



E^ifi) = El, - E, 



rbulk 
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where E^, is the total energy of the defect system, £', 
the total energy of the unperturbed system, Ua the num- 
ber of extra atoms of species a needed to create the 
defect, and rja the corresponding atomic chemical po- 
tential. Charge transition levels correspond to specific 
values of the electron chemical potential for which two 
charge states have equal formation energies. We consid- 
ered both thermodynamic and vertical charge transition 
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FIG. 2: (Color online) Comparison between charge transi- 



tion levels calculated with the semilocal (Jf' 



'^) and hybrid 



{Ji}^^^) functionals for a variety of defects in Si, SiC, Hf02, and 
Si02. The energy levels corresponding to the valence band 
maximum (VBM) and conduction band minimum (CBM) are 
also shown (squares). All energies are referred to a common 
reference level (j) (see text), shifted to coincide with the VBM 
in the hybrid scheme for convenience. For each material, A is 
the r.m.s. error with respect to the ideal alignment (dashed). 



levels. 

To compare defect levels in semilocal and hybrid 
density-functional schemes, it is necessary to use a com- 
mon reference level (f) external to the electronic sys- 
tem, i.e. defined on the basis of the nuclear potentials: 
71 = jj, + Ev — (f). This alignment scheme is illustrated 
in Fig. [ij In our formulation, we used the same pseu- 
dopotentials in the semilocal and hybrid calculations and 
trivially achieved such an alignment by taking cf) as the 
cell average of the local potential originating from the 
ionic pseudopotcntials. However, we note that our for- 
mulation does not imply any loss of generality and that 
a proper alignment can also be enforced when different 
pseudopotcntials are used. 

We calculated charge transition levels for the selected 
set of defects in Si, SiC, Hf02, and Si02 within both the 
semilocal and hybrid schemes. Charge transition levels 
obtained in either scheme were then aligned with respect 
to the common reference level (j) and reported in Fig. [51 
For each material, our results show that the defect lev- 
els calculated in the semilocal and hybrid schemes differ 
on average by at most 0.2 eV^when aligned in this way, 
despite the significantly larger concomitant variations ob- 
served for the band gaps. Since average shifts are sim- 
ilar in the four cases studied, the identified alignment 
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levels in both schemes 18 1. 

In the Slater approximation, the difference in charge 
transition levels can then be expressed as 
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FIG. 3: Optimal slopes derived from linear regressions of the 
data in Fig. [2] as a function of (a) experimental band gap and 
(b) average spread of the defect wave functions. 



is more impressive for large band gap materials where 
these shifts are small with respect to the band gap. In- 
deed, the average relative shift is only 2% for Si02, but 
increases to about 17% for Si. In particular, these results 
indicate that differences between charge transition levels 
are already well described at the semilocal level, lending 
support to alignment schemes in which the defect levels 
are anchored to experimental marker levels given by well 
characterized defects [J, Il7|. In this respect, a key result 
of our work is that the defect levels are positioned on an 
energy scale unaffected by band gap renormalization. 

To reveal systematic deviations with respect to the 
ideal alignment, we carried out linear regressions of the 
available data deriving optimal slopes [Fig. [3fa)] . Ideal 
alignment corresponds to unitary slope and is best il- 
lustrated for Si02 (slope of 1.08). When the band gap 
decreases, the optimal slope is found to increase indicat- 
ing that defects in the upper part of the band gap tend 
to follow the conduction band, while defects in the lower 
part of the band gap tend to follow the valence band. 
This tendency is most pronounced for Si (slope of 1.4). 

To provide a rationale for the obtained results, we first 
discard effects which affect the defect levels in a minor 
way. In all the cases studied, irrespective of band gap, 
the electron wave functions were found to be very simi- 
lar in the semilocal and hybrid schemes. Their effect on 
defect levels can be quantified by calculating total en- 
ergies at the hybrid level using electron wave functions 
optimized in the semilocal scheme without alowing for 
electron relaxation 



18[. In this way, we inferred that 



the differences between the defect levels calculated in the 
two schemes do not arise from variations of the electron 
wave functions. Since the observed trends also hold for 
the subset of vertical transitions, we discard effects as- 
sociated to differences in structural relaxation energies. 
The Slater transition-state approximation then allows us 
to focus on the single defect eigenstate rather than on 
the full manifold of occupied states [19|. This approxi- 
mation expresses the relevant total-energy differences by 
the energy eigenvalue of the defect state at half filling, 
/I = {ipY)\T-l\ip-D) — </^, and gives accurate charge transition 
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where the exchange potential V^ is evaluated at half fill- 
ing. Two different contributions can conceptually be 
distinguished in Eq. ([2]) referring to defect-defect and 
defect-bulk interactions. The separation is trivial for the 
Hartree-Fock exchange term, but requires some prescrip- 
tion for semilocal exchange. We first focus on defect- 
defect contributions. Since the defect wave function ipD 
is atomically localized, differences due to these interac- 
tions between charge transitions levels derived in the two 
schemes should be analogous to the corresponding differ- 
ences for ionization potentials and electron affinities of 
atoms and molecules. The latter quantities can be ex- 
pressed as total-energy differences and are already well 
described in semilocal approximations [20|,|2l|, as demon- 
strated by extensive quantum chemistry calculations [9|. 
Hence, this contribution is expected to give energy-level 
differences independent of the location of the defect level 
in the band gap. As for the defect-bulk contributions, 
it can be shown that they vanish in the limit of point- 
like defect states. When the defect wave function has a 
finite extent, these contributions depend on the degree 
of valence-band vs. conduction-band character of the de- 
fect state and can lead to a slope larger than 1 in Fig. [51 
To support this picture, we calculated average spreads of 
the defect wave functions in each material. Figure [3{b) 
clearly shows that the slopes of the linear regressions in- 
crease with these spreads, as the band gap decreases. 

Our results reveal a general trend which appears 
amenable to generalization. When the theoretical de- 
scription is improved, band edges in these materials 
undergo significant shifts but charge transition levels 
of atomically localized defects remain practically unaf- 
fected. This leads us to propose the following guideline to 
locate charge transition levels in the experimental band 
gap. First, ordinary semilocal density- functional calcu- 
lations are performed and charge transition levels of the 
targeted defect as well as band edges are determined. 
Second, the positions of the band edges are corrected 
through the use of a high-level electronic-structure the- 
ory which yields a band gap in agreement with experi- 
ment, e.g. through hybrid density-functional or GW cal- 
culations. This only requires a calculation for the bulk 
material, which is computationally less demanding than 
a defect calculation. Third, the defect level is located 
in the new band structure following the alignment pro- 
posed in this work. We demonstrate the applicability of 
this scheme for well-characterized defect levels: the two 
donor levels of substitutional Te in silicon fTem ) 12211 . the 
acceptor level of interstitial C in silicon (Ci) [23|, and 
the optical transition between the valence band and the 
E[ defect state in a-quartz 2J, |25| . The levels are first 



TABLE II: Comparison between calculated (/i"^) and mea- 
sured (/^°^P*) defect levels given with respect to the valence 
band. Theoretical levels are obtained from semilocal levels 
j-^scmiioc-j ^^JJ]-Q^gJ^ application of GW corrections to the va- 
lence band edge (Ae^^) [23]. Experimental data for the 
defects Tesi, d, and E[ are from Refs. [H, [H, and [H, 
respectively. 
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Defect 



q/q' 



Ae 



GW 



A* 



M 



GXpt 



Si Tesi 

Si Tesi 

Si Ci 

Si02 E[ 



0/+ 

0/- 

+/0 



0.2 
0.4 
0.5 
4.1 



-0.4 
-0.4 
-0.4 
-1.9 



0.6 
0.8 
0.9 
6.0 



0.6 
1.0 
1.0 
6.0 



determined with respect to the valence band maximum 
within the semilocal scheme 



26[. The position of the 



valence band is then adjusted according to recent GW 
calculations [27|. The resulting defect levels agree with 
the measured ones within the errors expected from our 
analysis f Table [IT)) . 

This guideline applies to atomically localized defect 
states and is clearly inappropriate for effective-mass-like 
defect levels which are tied to band edges. Application 
of this procedure also requires that the defect is well de- 
scribed already within the semilocal density-functional 
scheme. For instance, the defect level should fall within 
the reduced band gap of the latter scheme to preserve 
its localized nature. An inaccurate description may also 
result from the occurrence of competition between defect 
states featuring different degrees of localization 28 1. 

Our findings relate to other studies of defect levels and 
band gaps. Indeed, the band gap can also change as a 
result of a physical process, such as quantum confine- 
ment. For quantum dots of varying size, it has been 
shown that ionization potentials of deep defects remain 
constant as the band gap changes [29|. These poten- 
tials correspond to charge transition levels referred to 
the vacuum level. Another way to modify the band gap 
is achieved by changing the host material. It has been 
found that energy levels of transition-metal impurities 
within a set of isovalent semiconductors are aligned when 
referred to the vacuum level [30|. From the perspective 
of the present work, such an alignment is understood to 
the extent that the local chemistry of the defect is pre- 
served and a common reference potential can be identi- 
fied. Such transition-metal markers can then be used to 



predict band-offsets at interfaces J3l| . 

In conclusion, calculated energy levels of atomically 
localized defects generally remain tied to a suitably de- 
fined reference level as the description of the band gap is 
improved. This leads to a guideline for comparing calcu- 
lated and measured defect levels even when the adopted 
theoretical scheme is subject to the band gap problem. 
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